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Abstract 

The polaron has been of interest in condensed matter theory and field theory for 
about half a century, especially the limit of large coupling constant, a. It was not 
until 1983, however, that a proof of the asymptotic formula for the ground state energy 
was finally given by using difficult arguments involving the large deviation theory of 
path integrals. Here we derive the same asymptotic result, Eq ~ — Ca^, and with 
explicit error bounds, by simple, rigorous methods applied directly to the Hamiltonian. 
Our method is easily generalizable to other settings, e.g., the excitonic and magnetic 
polarons. 

The polaron Hamiltonian of Frohlich is a model for the Coulomb interaction of one or 
more electrons with the quantized phonons of an ionic crystal. In the course of time it was 
also seen to be an interesting field theory model of non-relativistic particles interacting with 
a scalar boson field, and it was widely studied p|, H, ^, H, ^, ||, 3 in both contexts. 



©1996 by the authors. Reproduction of this article, by any means, is permitted for non-commercial 
purposes. 



1 



The model has one dimensionless couphng constant, a, and it was noticed very early that 
there seems to be a qualitative difference in the ground state between the wealc coupling 
regime, well described by perturbation theory in which the electron is spread out, andthe 
strong coupling regime in which the electron appears to be trapped in a phonon hole of its 



own making. (This question of trapping seems to have first been considered by Landau |]T0 
in the context of a classical phonon field.) 

By now it seems doubtful that such a trapping actually occurs |ll|, |T2|, but it remains 
true that the calculation of the ground state energy, £'o(a), is very different in the two 
regimes. The strong coupling theory was studied by Pekar (see also H, il i|, |^) who 
hypothesized that in this limit the total ground state wave function could be taken to be 
a product of an electronic function 0(x) and a phonon function \^). It is a fact that this 
ansatz, called the adiabatic approximation, yields exactly the same result for the ground 
state energy as the corresponding model with classical phonons. This stems from the fact 
that the phonon operators enter only linearly and quadratically in the Hamiltonian. 

Frohlich's Hamiltonian in appropriate units is 



— e + h.c. 
k 



(1) 



where p = — iV is the electron momentum operator, x is its coordinate, and V is the 
volume of the crystal, which tends to infinity. The k's are the usual normal modes ( e.g., 
k = 27rV^~^/^(ni, n2, ns) for a cubic box ) and, as usual, V'^Y^k (2vr)~^/ (Pk. Also, as 
usual, the k = mode is omitted. When (\E'|if|\E') is computed with Pekar 's ansatz one easily 
finds that the phonon part |^) can be easily evaluated by 'completing the square' and the 
resulting minimizes the energy 



^(0) = J \V(j)\'d'x -a J J ^^^^^^d'xd^y (2) 

subject to / (f){'x)^d^x = 1. This, as we mentioned above, is what one would get if the Ok's 
are replaced by c-numbers ^k- 

It has been proved |T^ that there is exactly one 0, up to translations, that minimizes 



S{(j)). By scaling, the minimum energy (call it Ep{a)) is proportional to — a^. The question 
we address is this: Is this energy, Ep{a), asymptotically exact as a tends to infinity? Since 
it is given by a variational calculation, Ep{a) is an upper bound to Eo{a); the problem is to 
find a lower bound that agrees with Ep{a) to leading order. We provide such a lower bound 



in ( p7| j below. 

For a long time more or less the only thing that could be done to validate Ep{a) rigorously 
was a lower bound |0 of the right order, — a^, but which was about a factor of 3 too large. 
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It was not until 1983, with a precursor in |T^, that a complete proof of the asymptotic 
correctness of Ep{a) was given by Donsker and Varadhan The proof starts with a 

familiar expression for the ground state energy, 

i?o(«)= lim ^ln(vl>|e-^^|vl>), (3) 

with \l/ such that its spectral resolution contains the ground state energy or low energy 
spectrum of H but is otherwise arbitrary. Now in fact there is a relatively simple Feynman- 
Kac formula for the semigroup e~^^ first noted by Feynman in [Q; by integrating out the 
phonon coordinates in this formula, one finds that the right side of (j^) can be replaced by 
the expression 

1 f /-/^ /■* e"(*"^) 1 
lim — InE'Kexpfa / / — -—-dsdt)>. (4) 

/3-oo/? \ io Jo |x(t)-x(s)| '] ^ ' 

Here, E represents path space expectation for the Brownian motion x(t). The Feynman-Kac 
formula and functional integral are discussed in some detail by Ginibre [|I^ and Roepstorff 
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It is this limit that Donsker and Varadhan analyze, in effect showing that the Pekar 
functional (Q) provides a kind of large deviation rate function for the functional integral in (Q). 
Their proof is very far from simple and requires an enormous knowledge of large deviation 
theory of Brownian motion, a subject unknown to most physicists, and it does not generalize 
easily to other models. Owing to the double integral in (H) it is necessary to construct a large 
deviation theory for "infinitely many variables". Moreover, no error bounds were provided 



17| . We note that previously Feynman had used this formula to estimate the ground state 
energy from above, but his estimate is not sharp even for a — >■ oo. See also where this 
functional integral is again employed to obtain a lower bound on the ground state energy 
which, however, is not sharp. 

Regarding the matter of error bounds, we note that Gross takes advantage of momen- 
tum conservation to eliminate three phonon degrees of freedom and transforms the Hamilto- 
nian H to a new form involving a sum of a^-, a*^- and a~^-terms. (He continues to employ 
the Pekar ansatz.) This form of the Hamiltonian perhaps suggests that the ground state 
energy itself has an asymptotic expansion in inverse powers of a^. But the terms in this 
transformed Hamiltonian are unbounded operators, and so the order of their contributions 
to the energy is by no means obvious. The bound we obtain on the error is actually 0{a^^^), 
which is to be compared with the main term O(a^). 

In any event, there is the following problem: The Pekar ansatz is based on the physi- 
cally appealing notion that at large coupling the phonons cannot follow the rapidly moving 
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electron (as they do at weak coupling) and so resign themselves to interacting only with 
the 'mean' electronic density, 0(x)^. What, exactly, do these somewhat anthropomorphic 
words mean? If they are so very physical, they should be quantifiable and it should not take 
four decades to build a proof of their correctness. Moreover, one should be able to find a 
proof that is relatively simple (since the idea is a simple one), and one that also yields some 
kind of quantitative error estimate. It should also be robust enough to allow an easy exten- 
sion to some variations of H, such as having several electrons instead of one (the polaronic 
exciton), inclusion of magnetic fields, an electron kinetic energy, T(p), other than (for 
which the semigroup exp[— T(p)] might not even have a positive kernel), accommodation 
of k-dependence and nonlinearity in the phonon self energies, electron-phonon interaction 
energies other than l/|k|, etc. 

The method presented here satisfies, we believe, the criteria of simplicity, robustness and, 
of course, rigor. For this reason it might be of general interest in condensed matter theory or 
field theory. Despite its generality we shall, for clarity, restrict our discussion to the original 
Hamiltonian, H . 

Before going into the details, it is helpful to give an overview of our method. The first 
thing to notice is that Pekar's ansatz amounts, mathematically, to saying that we can replace 
the operators Ck by c-numbers ^k- One might think, at first, that such a replacement always 
leads to a lower bound for the energy, but this is false. (If it were true, the small a energy 
would have to be proportional to —a^ instead of the much more negative, correct value —a.) 
To pursue this idea, nevertheless, the natural tool to think of is coherent states for each 
phonon mode k. As is well known, these states are indexed by complex numbers ^ = p + iq 
and, if \^){^\ is the projector onto the coherent state |^), we have that 

/ \0{^\d^dC^I, (5) 

J^\0{^\d^dC^a (6) 

and 

Jm'-m){C\dCdC^a*a. (7) 

It is the eoctra term —1 in the last integral that kills the lower bound and gives an unwanted 
energy —1 for each phonon mode (and hence a negative contribution to the energy that is 
infinite both because there is no infrared and no ultraviolet cutoff). In other words, coherent 
states would give us what we want (effectively replacing the operator a by the number ^), 
were it not for the unfortunate fact that the positive operator a* a is represented by the 
nonpositive symbol, — 1. 
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Our remedy will be to reduce the effective number of plionon modes to a finite number 
of 0(a^/^), independent of V. This will not sacrifice rigor because we shall prove that the 
reduction affects the energy only to 0{a^^^) at most. These modes will be quite different 
from the original modes; indeed they will be the Ok modes summed over boxes in k-space 
of size a^/^. Thus, our physical description of strong coupling will be a little different from 
the conventional one. Instead of saying that the phonons cannot follow the electron, our point 
of view will be that the electron significantly excites only finitely many field modes. 

This mode reduction is accomplished in several steps which can be outlined as follows. 

• I. Using a simple commutator estimate (as in Schrodinger's elementary proof of Heisen- 
berg's uncertainty principle |Q) we can show that |k| values larger than K ^ a^^^ can be 
ignored with an energy cost of only a^^^. This eliminates the ultraviolet problem, i.e., the 
fact that |k| is not summable. 

• II. With the same energy error we can localize the electron to a cube of side length 

• III. We decompose the ball |k| < K into blocks of size a^/^, each containing n k-values 
(with n ~ Va^/^). There are about a^^^ blocks. Because the electron coordinate x has been 
localized, the function exp(ik ■ x) in if can be replaced, in a block B, by exp(zkB ■ x), where 
ks is any conveniently chosen point in B ; The energy penalty is again a^^^. The effect of 
this decoupling is that the electron interacts only with the one mode Ab ~ ""-^^^^ Z^keB '^k, 
from block B. Thus, the electron ends up interacting with only a^^^ modes , which is a finite 
number! 

• IV. These active modes can now be represented by coherent state integrals, as above. 
The unwanted —1 term now contributes only a^^^ to Eq, that being the number of modes. 

The details will now be given. 



I. We first establish a commutator inequality which will show that the large k modes of 
the Hamiltonian may be discarded at the price of only a small decrease in the coefficient 
of the term of the Hamiltonian. The argument given here is incorrect. Please see the 
erratum at the end of the article. For any normalized state with expectation (■), we have 

|fc.ll(«ke^'-'')l = l([P.,«ke''-])| 

< 2(p|)V2(a*aO^/^ (8) 

for each phonon momentum coordinate kj, j = 1, 2, 3, just by the Schwarz inequality. Squar- 
ing this inequality, summing over j, and then taking the square root, we obtain 

l(«ke^'-^)l<^(P^)^/'Kak)^/^. (9) 
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A similar inequality holds for Ok replaced by a^. It follows that for any e > 0, and with 
K = Sa/ne, 



,47ra 
' V ' 



E 

k|>iC L 



|k|>iC 



|k|>_ft: |k|>-K: 

|k|>X 

< £(p2) + ^ ( a*ak ), 

|k|>X 



(10) 



again by the Schwarz inequality. (We have taken the limit — >■ oo in the sum.) 

The above inequality ([10|) is an ultraviolet bound. It shows that the Hamiltonian H is 
bounded below by a new one Hk-, i-e., H > H^, with 



H 



K 



^1 ^)P + L «k«k + (-r^)^ L 

|k|<i<' ^ |k|<A' L 



TlK' 



— e + h.c. 
k 



By completing the square, we have that 

,47ra 1/2 



V 



^e^*^-^ + h.c. 
k 



> 



Ana 



:i2) 



and if K is chosen so that the coefficient of vanishes in Hk-, i.e., K = ^a/Ti^ then 

16 o 



*^ |k|<K 



^« 



(13) 



showing that H is indeed bounded below by O(a^). But to obtain the sharp lower bound, we 
will take a much larger iiT, namely K = {8/nci)a^^^. Here and below, the q's are constants 
independent of a. 



II. The next step is to localize the electron. Although a-priori the electron is in the box of 
volume V, it actually can be confined to a box of much smaller size, with only a slight relative 
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increase in its energy. More precisely, we have the following: Let \E' be any normalized state 
of the electron and phonons and define E = {"^IHkI"^)- Then, given AE > 0, there exists a 
function of the electron coordinate x alone (but depending on with support in some 
cube of side length L = n{3/ AE)^^"^ such that 

{(j)<if\HK\(j)'i/) / {(jy^lcjy'if) < E + AE. (14) 

To see this, let 0(x) = n^=i cos((AE/3)^/2a;j) inside the cube of side length 7r(3/AE)^/2 
centered at the origin and = outside the cube. Let 0y(x) = 0(x — y) and consider the 
integral 

= |(|V0p-AE|0p)rf3^ = O. (15) 

(Note here that the cross terms V0yV\E' vanish after the ^/-integration.) Evidently, there 
must be a point y such that at this point, the integrand on the left side of (|r^) is non- 
positive (and where (0y\l/|0y^) is non-zero). This (py is the (p that we need. The electron 
localization is now complete. 

For our purposes, we take AE = C2a^^^, and the assertion above together with the 
ultraviolet bound then implies that the ground state energy Eq of H satisfies 

Eo > inf - caa^/^ (16) 

where the infimum is taken over all normalized \l/'s but having their x-support in a cube of 
side length no larger than L = (3/c2)^^^vrQ;~^/^° somewhere in the large volume V. 

III. The next step is to group the phonon modes together into blocks, which we take 
to be cubes with sides of length P = c^a^^^ or, more precisely, the portion of those cubes 
lying in the big ball {k : |k| < K}. Let 5 be a block (cube) of momenta with sides of 
length P = csa^/^ and let be any fixed point within this block whose precise position will 
be determined later. It is clear that the set of momenta {k : |k| < K} can be covered by 
N = ^nK^/P^ + lower order = (2^^/(37r^cfc|))a^/^ -|- lower order such blocks. 

Assume that x itself varies over a cube of side length no bigger than L with center 
at a point xq which we take to be the origin 0. (If xq is not the origin, phase factors 
exp (±i(k — k^) ■ Xq) standing before Ok and in the Hamiltonian Hk are readily eliminated 
by a unitary transformation involving the phonon variables only.) We have that for k G 5, 

Jg^k.x _ g.k^.x| < |(k - ks) ■ X| < ^(i)V2^-3/10_ ^17^ 

2 Co 
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Then, for any 5 > 0, we have the following inequality, which is obtained by completing the 
square, as in (|T^) . 



EE 

B keB 



- 5V \ 2 ^C2^ ; I I 



|k|<_fs: 



CiC2(5 



We make the choice 6 = Ci^a~^l^ . 

The above inequality (^) permits us to bound Hk (still restricted to states with x-support 
in a cube of sides with size V) from below in the following way: 



Hk > (l-Cia-i/V + EE[(l-C4«" 



l/5^ 



B keB 



+ 



V 



108c|«9/5. 



C1C2C4 



(19) 



Note that, within a single block, the exponential factors multiplying the at- and a^-terms 
are now independent of k. With this in mind, we define, for each block B, a block annihilation 
operator 

Ok 



A 



B 



E \^n-'^' E 

keB 1^1 



(20) 



keB 



and its corresponding adjoint A*^. These operators are properly normalized boson modes, 
i.e., 

[Ab,A*s,]=5b,b'. (21) 

Then 

E <«k > A%A^, (22) 

kS-B 

and so the operator terms on the right hand side of ([T9|) exceed the operator H'^"'^^ = 
^Wocfc^|]^^|) defined by 



H 



block 
K 



:i-cia-vv+E 



[l-C4a-^^^)A*sA 



B 



+ (-^E |k|-')^/2(ABe^'^-''+/i.c.) 



keB 



(23) 
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Referring to (|T6|) and (|19|), we can summarize the situation: The ground state energy Eq of 
H satisfies 

Eo > inf sup {^M^^'^m - (C2 + ^^)a^/^ (24) 

* {kg} C1C2C4 

where the infimum can now even be taken over \I/'s with no restriction on their x-support, 
and with the k^'s chosen optimally, depending on 

IV. In order to get a lower bound on H'^°'^^, we use the technology of coherent states. Let 
1^) = vr~^/^exp (— + ^A*)|0) denote a normalized coherent state for a single harmonic 
oscillator. In terms of these states one has the identity 

exp(M) exp{uA*) = J exp (/.^ + uCm{^\d^dC, (25) 

from which the identities for the operators A and A* A mentioned in the introduction are 
readily obtained. 

Let \l/ be any normalized state and let |^) = \^b) denote a tensor product of coherent 
states of the block modes corresponding to the operators {As}- Set ^^^(x) = {^\'^)phonon, 
where the inner product, ( | )phonon, is just over the phonon variables, not x. Then 

(Xl/|ij^-'=|vl/) = J {^^\h^\^f^),i,,,ronY[Bd^BdeB, (26) 

where the inner product in the integrand, ( | ) electron, is over the electronic coordinates x, 
and is the Schrodinger operator 



;i-cia-^/V + E 

B L 



+(^Elkr^)^/^(eBe^''-^ + c.c.) 



1- ca-'/'mnl' - 1) 

(27) 



Now the supremum over of (|26|) exceeds what is obtained by completing the square, 
namely, 

(28) 

where p^(k)= (^^| exp (ik ■ ^)\'^^) electron is the Fourier transform of |\&^p(x). (Here, we can 
remove the |k| < restriction on the k-sum.) We emphasize again that were it not for 
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the grouping of the phonons into blocks, then this coherent state estimate would contain a 
negative constant equal to the number of phonon modes with |k| < ii', which is infinite in 
the V ^ oo limit, instead of = 0{a^^^). 

At this point we choose to be a point in B where the function of k given by 

/^n.«. (29) 

is minimal. This point depends on \E'. With this choice, (^) is, in turn, larger than 

/ [(1- C,a-'/')mp^\^,W,ron - (i_^'j-V5)y E J^^] H.^B-^^- (30) 

But the integrand of this integral is {'^ ^\'^ ^) electron x S{-^^/{{-^^\-^^)ll^^t^^J) , where S{-) 
is just the Pekar functional (^, (but with coefficients altered by 0(a;^^/^)); the k-sum is 
the Fourier series for the Coulomb self-energy term of S. Thus, the integral (^) exceeds 
(1 — Cia~^/^)~^(1 — C4a~^/^)~^i?p(a), with Ep{a) being the Pekar minimum energy of (0). 

Taking into account the altered coefficients in the Pekar functional, the lower bound ([2^), 
and the definition of N, we obtain 

inSr^ oil 

Eoia) > Ep{a) - (c2 + ^^ + -^^ + c^)o?l'> -o{o?l'>) (31) 

C1C2C4 dTT^CfCs 



where Ep{a) = —cpo?^ cp = 0.109,[jTTl, and where C5 = (ci + 2c4)cp accounts for the altered 
coefficients in the variational principle. This completes the proof of the lower bound. Opti- 
mizing the coefficients, we find that the a^^^-error term in (^) is no greater than 3.822a9/^ 
This work was partially supported by NSF grant PHY95-13072(E.H.L.). 
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Erratum 

We are grateful to Professor Audrey V. Soldatov of the Moscow Steklov Mathematical 
Institute for calling our attention to an error in our paper. The commutator inequality (8) 
in our step I, namely |/cj||(ake*'*"'')| < 2(p^)^/^(akak)^/^, is not correct. Rather, the right 
side of this inequality should be {p]y^'^{{at(^ky^^ + (Ok'^k)^^^) or a related expression. The 
extra factor (Ok'^k)^^^ with the a^^ and aj^ not in normal order generates uncontrolled mischief 
with, for example, the right side of the ultraviolet bound (10) containing an additional term 
E|k|>i^ 1/2 = oo. 

The situation is remedied with the help of the method introduced by Lieb and Yamazaki 
[14] to obtain the previous rigorous lower bound on the polaron energy. Our main result, 
(31), is still valid. Indeed, it is improved slightly. 

Define the (vector) operator Z = {Zi, Z2, Z3) with components 

= (^)^/^ E O"^' ^' = ^'2' 3- (32) 
^ |k|>ii: 1*^1 
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Then the commutator estimate (8) is replaced by 



V 



\k\>K 



< 2(p2)V2(_(Z-Z*)2)V2 < 2(p2)V2(2(Z*Z + ZZ*))^/2 

< £(p2) + -(Z*Z + ZZ*). 



(33) 



Now, each component Zj can be thought of as a single (unnormalized) oscillator mode having 
commutator with its adjoint, [Zj,Z*] = (47rQ;/l/) X]|k|>K — > 2a/3nK; moreover, Z^ 
and Z* commute for i ^ j (i.e., these modes are orthogonal). Using these facts, we have that 



2 4 2 / 20- 

-(Z*Z + ZZ*) = -(Z*Z) + -( — 
e EE K'kK 

< E («k«k) + 3/2 

\k\>K 



(34) 



if we choose e — 8a/37rK, which is smaller and better by a factor 1/3 from the e in the article. 
Here we have employed an orthogonal rotation of coordinates bringing I]|k|>x '^k'^k i^^o ^ 
form (4/£)Z*Z+non-negative operators. (Compare Eqs. (21,22) of the article.) Combining 
these inequalities, we obtain 



,47rQ; 



E 

ik|>i^ 



^e*-'') + c.c. 
' k 



E< 



) + 3/2. 



(35) 



This last inequahty is our replacement for the ultraviolet bound (10). It follows that 
H > Hk — 3/2 where is as in Eq.(ll), but with the coefficient of p^ given by (l — Sa/SnK) 
rather than (1 — 8a /nK). With the choice K = 8a/37r, inequality (13) becomes H > 
— (16a^/37r^) — 3/2, a bound at least consistent with a known upper bound for the ground 
state energy linear in a. 

The remainder of the article is an analysis of and needs only minor modification. 
The coefficient of p^ in Eqs.(19,23,27,28,30) should be (1 - and, at the end of the 

article, C5 = (ci/3 + 2c4)cp. Due to the smaller value of e defined above, our estimate on 
the coefficient of a^^^ in (31) is slightly improved to 2.337, rather than 3.822 as reported. 
Of course, our lower bound for the ground state energy is decreased merely by the constant 
—3/2, which is unimportant on a scale of a^^^. 
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